use 	"./output/cities_expandedpatents", clear

keep if f_myr == tm(1918m9)

** PLOT OPTIONS **
set scheme plotplainblind

histogram days_npis, w(5) freq start(25) ///
	xlabel(, nogrid) ylabel(, nogrid) fc(gs15) ///
	xtitle("Days of NPIs (5-day bins)") ytitle("Count of Cities") 
	
graph export "$RES/FigureAX_histogramNPI.png", replace